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Abstract 

A simple model is proposed for the buckling and coiling instability of a viscous "fluid rope" 
falling on a plane. By regarding a fluid rope as a one-dimensional flow, this model accounts for 
only the axial and shared viscous forces. Our model successfully reproduces several experiments 
with no adjustable parameters, such as the existence of three distinct coiling regimes reported in 
Phys. Rev. Lett. 93, 214502 (2004). Our model allows for the discussion of unsteady motion. 
An expression for the critical fall height at which the coiling frequency changes from a decrease 
to increase was phenomenologically derived. It was found that the coil-uncoil transition shows 
remarkable hysteresis only for weak gravity condition. 

PACS numbers: 47.20.Gv 
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In the low Reynolds number regime, the dynamics of viscous fluid confined in rigid 
boundaries do not exhibit complex behavior because the possible eigenmodes are limited. 
However, if the fluid interface can move freely and deform largely, this is not the case, and 
rich variety of dynamics can emerge below a critical Reynolds number. Backing and coiling 



instabilities of fluid jets are exam 
decades in the laboratory 



re examples of ; 

Q,y,y,y,y,i6]. 



such systems, and have been studied for several 
A fluid rope and coiling can also be observed in 
daily life, such as when honey is poured from a teaspoon onto toast. 

Although the mathematical treatment of a largely deformed fluid is not theoretically 
straightforward, the conditions for the onset of coiling instability have been thoroughly 
examined in terms of linear stability analysis with some simplifications and assumptions 

m. 

Recently, Ribe derived the differential equations to describe steady state coiling for a 



very thin rotating rope, and showed the numerical solutions for various conditions [10|. His 
analysis reveals that the steady solution is multi-valued so that there might be a discontinuity 
in the selected coiling frequency. 

Ribe also demonstrated the existence of three distinct coiling modes; viscous, gravita- 
tional, and inertial regimes. Phenomenologically, one could understand that buckling and 
coiling instability occurs under the mechanical balance between the driving force of a steady 
flow and the internal stress due to viscosity. Signifying the magnitudes of force for fluid 
injection as Fp, the gravitational force acting to fluid as Fq, the inertial force as Fj, and the 
yielding stress due to viscosity as Fy, the viscous coiling regime corresponds the condition 
Fp ~ Fy, the gravitational coiling regime Fq ~ Fy, and the inertial coiling regime Fi ~ Fy, 
respectively. Dimensional analysis allows for the estimation of the coiling frequencies for each 
regime as 

where Q is the flow rate, H is the fall height, v is the kinetic viscosity, g is the gravitational 
acceleration, and a is the radius of the rope. 

Although these studies seem to succeed in giving the onset of coiling with infinitesimal 




FIG. 1: A schematic view of a fluid rope coiling. 



amplitude and steady coiling frequencies, in order to discuss the stability and further turbu- 
lent states, a numerical model is still needed that can describe the entire dynamics including 
transient states with fewer computational costs. 

In this Letter, an alternative numerical model for a small deflection in a fluid rope is 
proposed and coiling frequency as a function of fall height is investigated. In addition, the 
coil-uncoil transition and its hysteresis effect are discussed, as well as the transition height 
from viscous to gravitational coiling. 

Figured] is a schematic view of a fluid rope injected from an orifice at a sufficient height. 
We fix the origin of the reference flame at the point where the rope begins to coil and align 
the z axis parallel to gravity g. In the present model, we require a "steady coiling state" 
in which the upward growth speed of the coil is balanced by the speed of downward flow at 
the top of coil, i.e. the point at which the rope changes to the coil does not move. 

We consider a uniform flow within the rope and write the flow velocity as w = wn, where 
n = (n x , n y , n z ) is the unit vector tangential to the rope. Since we restrict the present model 
to the case of slight deflection, the vector n is almost parallel to the z axis. We take account 
of only the first order of n x and n y , and replace the derivative with respect to n by the 

stress is written as a n = 3r]dw/dz, 



ll|. Hereafter, we scale flow speed 



derivative with respect to z. In this manner, the axia 
where 3rj is the extension (or compression) viscosity 
with the injection speed at the nozzle W[ Q and length with the orifice diameter d. Then, the 
dimensionless equation of motion for w is 

8 8\ 3d f a dw\ 1 



dt dz J SRe dz \ dz J Fr 



where Re = d\w\ n \/u is the Reynolds number and Fr = wf n / gd is the Froude number. S is 
the rope's cross-sectional area taken to be parallel to xy plane. We determine S by solving 
the following continuity equation 

f = !<*■>• < 5 > 

In order to describe the deformation of the fluid rope, we divide the rope into thin cylindrical 
elements whose motion is restricted in the xy plane, as shown in Fig. [TJ Denoting the 
position of the element q = (q x ,q y ) and velocity u = (u x ,u y ), these quantities obey 

+ ™7jI ( 6 ) 



dt dz J 

where i = x and y. The total stress acts on the cylindrical element is the sum of viscous 
shear stress Oi = i]d(ui + wni)/dz and the iih component of the axial stress Ui<j n . Therefore, 
we obtain the equation of motion for a cylindrical element as 

d d\ 1 d [ ( dui dw drii \ ] . . 

m +W Y Z ) Ul= SRed- Z { S {-dz- + An ^ + W ^z-)l (?) 

Next, we discuss the boundary conditions of the present model. At the injection point 
z = H, the fluid rope is fixed, or 

(&(#) = 0, Ui(H) = 0, (8) 

and S(H) — 7r/4. We neglect the relaxation of Poiseuille flow to plug flow in the neighbor- 
hood of the orifice, and assume a uniform flow at the injection point as w(H) = —1. Let us 
turn to the boundary conditions at z = 0. We assume that 

q'M = 0, <(0) = 0, (9) 

where a prime indicates a derivative with respect to z. Let w max be the maximum of 
Experimentally, Cruickshank and Munson found a relation Wb ~ — 0.13w max , where u>& is 
the downward flow speed at the top of coiling region [3] . This relation is fairly independent 
of flow conditions such as the viscosity of the liquid, the orifice diameter or the fall height. 
The velocity w(0) in the present model corresponds to Wj,, we therefore require the following 
time dependent boundary condition 

iw(0) = -0.13uw (10) 
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FIG. 2: Trajectories of the bottom of the model rope starts at t = 0: (a) Trajectory in the case 
that axial flow is stable at the steady state (Re = 10.0, Fr = 100.0, H/d = 10). (b) Trajectory of 
circular coiling (Re = 1.0, Fr = 1.0, H/d = 10). 



We numerically solve the equations (BJ, ([5]), (EJ) and (J7D setting an axial flow, w(z) = — 1, 
S(z) = 7r/4, qi(z) = 0, and Ui(z) = as an initial condition. For qi(z), a small roughness 
with amplitude of 0.01 is given. The tangent vector n is calculated from the derivative of 



Most of experimenta 
of the fall height [l, 2, 3 



q with respect to z as n = (dq x /dz, dq y /dz, 1) jC with C = ^(dq x /dz) 2 + (dq y /dz) 2 + 1. 
The control parameters of the numerical calculations are Re, Fr and the fall height H. 

Figure [2] shows trajectories of q in the plane z = 0. Simulation movies are available [ljj. 
In the case that Re is sufficiently large, the axial stagnation flow is stable, thus the trajectory 
converges to the origin, as shown in Fig. [2](a). Such an axial flow could appear under no 
gravity condition. When we set Re smaller than a critical value, the flow starts to oscillate 
and a circular trajectory would appear at the steady state (Fig. EJ^b)). 

works on fluid rope coiling discuss coiling frequency as a function 
, U, [5| or the rope's radius [6(. Let us compare the coiling frequency 
of the present model with the experiments. To investigate the viscous coiling regime as a 
first step, Fr is set to a very large value, since the viscous coiling takes place when the effect 
of gravity is small. Figure [3](a) shows a dimensionless frequency VLd/w in as a function of 
fall height. We found that Q oc H^ 1 for lower fall heights, which agrees with the scaling 
relation in viscous coiling. For H/d ~ 10, the frequency is almost constant. This behavior 
corresponds to the gravitational coiling, because the rope's radius hardly depends on the 
fall height. In figure [21(b), the frequency-height curves are rescaled with Qy and 111 
order to compare the results with the experiment by M. Maleki et. al. (see the inset). 
The transition from viscous to gravitational coiling occurs at Qg/^v — 1- These results 
are in good agreement with the experiments including the transient regimes. For a higher 
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FIG. 3: Transition from viscous to gravitational coiling for Fr = 100, 200 and 400. (a) Dimension- 
less frequency £ld/w- m versus height in condition of Re = 3. (b) Frequency-height curve rescaled 
using the fly and Oq. The inset shows the experimental result by M. Maleki et. al. 4|. 



fall height of H ~ 10 2 , the frequency increases as Q oc H, which agrees with the earliest 
observation by G. Barnes and J. MacKenze l|, |2[]. 

Supposing strong stretching a <C d due to gravity, the scaling relation Eq. (j3J) pre- 
dicts Q oc if 10 / 3 [5]. In this case, we reasonably choose a characteristic time scale y/d/g. 
According to simple dimensional arguments [3], frequency Qy/d/g must have the form: 



(ii) 



where v = rj/p is the kinematic viscosity. Here it is noticed that the first parameter in this 
function is the ratio of the viscous force to the gravitational force Fy/Fc, the second is the 
ratio of the gravitational force times the inertia force to the viscous force FqFi / Fy 2 , and the 
third is the slenderness ratio. Because the strong stretching condition requires Fy/Fc < 1 
and H/d ^> 1, the inertial coiling regime appears in the condition: 



uQ gd 3 ^ H 
" < 1, — > 1 



d 



(12) 



Experimentally, Cruickshank and Munson first found Q oc H 2 for H/d higher than 10, in 
the condition of vQ/gd 4 ~ 0.1-1.7, gd 3 / 
gd 3 /v 2 ~ 10" 3 and vQ/gd' ~ 0.5 



a. 



.1 3. 



Ribe also observed Vt oc H for 



As shown in Fig. HJ the present model gives Q oc if 2 2 for the similar conditions. These 
experimental and theoretical results are successfully consistent. On the other hand, in the 
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FIG. 4: Dimensionless frequency Q^d/g versus height obtained from the present model. Parame- 

ters are: vQ/gd 4 = 0.2 and gd 3 jv 2 = 0.2 (open circle), vQjgd 4 = 0.3 and gd 3 jv 2 = 0.3 (triangle) 
and uQ/gd 4 = 0.2 and gd 3 ' jv 2 = 0.5 (rectangle). The inset shows the experimental result by M. 



gravitational to inertial transitional range, the experiment shows a remarkable discontinuous 
jump in the frequency (inset in Fig. Hj). The present model, though, does not reproduce this 
behavior. 

Next, let us discuss critical fall height H*, around which coiling frequency changes from 
decrease to increase (transition from viscous coiling to gravitational coiling). In figure [3](a), 
the transition occurs at H* ~ Gd. However, the H* should generally be a function of viscosity 
and gravitational acceleration. For the first step, we focus on the position z = ( a,t which 
the flow velocity w reaches a maximum. Note that, internal stress of fluid rope changes 
from tensile to compressive at this point. After a flow reaches steady state, it is observed 
that ( hardly depends on time, even though the rope is oscillating. In Fig. [5](a), we plot ( 
at steady state as a function of fall height. For a shorter fall height, the ( is equal to H. 
This indicates that the maximum of w appears at the point of injection and the fluid rope is 
wholly compressed. Comparing ( with the frequency Vtd/w- m shown in Figj5](b), we can find 
that the frequency decreases only when the relation ( = H occurs. Therefore, viscous coiling 
appears in the case in which the whole the rope is compressed. Because this feature can be 
seen for a wide range of Re and Fr as long as the viscous coiling regime exists, we claim that 
the H* is identical to the maximum of (. The maximum of ( can be realized as a relaxation 
length of which the effect of boundary condition at z = can travel through a rope. Using 
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FIG. 5: (a) Numerically calculated value of z c for Re = 3 and Fr = 100. The solid line is the plot 
of £ = H . (b) Same data with Fig. [3] for comparison. 
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FIG. 6: (a) Numerically calculated coiling radius for Fr = 10 4 and H = 15 with increasing and 
decreasing Re. (b) The coiling regime in Re-Fr plane for increasing (A) and decreasing (y) 
Reynolds number. 



g and z/, we can uniquely construct a dimension of length as g 1 / 3 z/ 2 / 3 . Therefore, 

H* = max(C) oc g- 1/3 v 2/ '\ (13) 

which is fairly close to the numerically calculated exponents max (£) oc g-°- 36 i/°- S6 . 

Finally, we refer to a new and unforeseen result obtained from the present model. It 
is empirically known that fluid ropes do not coil if the Reynolds number of the flow is 
sufficiently large. We investigated the upper limit of Reynolds number of the present model 
and found that the transition occurs with remarkable hysteresis. In the simulation, we 
continuously add perturbation of amplitude 10 _6 <i for gj. Fixing Fr and H, the value of Re 
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is slowly increased or decreased not faster than 0.1% change per one cycle of coiling. Figure 
Efa) shows the coil-uncoil transition for weak gravity condition (Fr=10 4 ). For increasing Re, 
the coiling radius sharply falls down at Re ~ 3.3(= Re ; *) and for decreasing Re, the radius 
somewhat gradually rises at Re ~ 0.7(= RejQ. Axial flow is absolutely stable at Re > Rej*, 
and is absolutely unstable at Re < Re z *. In the intermediate region (Re z * < Re < Re£), 
fluid rope starts to coil if the amplitude of the perturbation is of order unity. In Fig. 0(b), 
we show the coiling regime in Re-Fr plane for increasing and decreasing Reynolds number. 
The hysteresis suddenly disappears for Fr less than about 400. We do not have a clear 
explanation for this drastic behavior yet. However we observe that the height ( starts to 
decrease from H roughly at Fr < 10 3 , therefore, a crossover of viscous and gravitational 
coiling regime should exist around Fr = 400. From these facts, we expect that the viscous 
coiling has hysteresis effect for coil-uncoil transition. To verify this, further experiments or 
direct numerical simulation should be carried out. 

The model we proposed in the present Letter could also reproduce recently reported 
meandering instability of viscous fluid rope falling onto a moving surface by simply 
replacing one of the boundary condition in Eq. ([9]) as q' x (0) = s, where s is the constant 
share rate. Further analysis of this problem is our next work and will appear elsewhere in 
near future. 

We would like to thank Dr. Katsuhiko Sato and Dr. Toshihiro Kawakatsu for their 
critical comments. 
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